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We outline the confluence of three novel theoretical fields in our modeling of Gamma- 
Ray Bursts (GRBs): 1) the ultrarelativistic regime of a shock front expanding with a 
Lorentz gamma factor ~ 300; 2) the quantum vacuum polarization process leading to an 
electron-positron plasma originating the shock front; and 3) the general relativistic pro¬ 
cess of energy extraction from a black hole originating the vacuum polarization process. 
There are two different classes of GRBs: the long GRBs and the short GRBs. We here 
address the issue of the long GRBs. The theoretical understanding of the long GRBs has 
led to the detailed description of their luminosities in fixed energy bands, of their spec¬ 
tral features and made also possible to probe the astrophysical scenario in which they 
originate. We are specially interested, in this report, to a subclass of long GRBs which 
appear to be accompanied by a supernova explosion. We are considering two specific 
examples: GRB980425/SN1998bw and GRB030329/SN2003dh. While these supernovae 
appear to have a standard energetics of 10 49 ergs, the GRBs are highly variable and can 
have energetics 10 4 10 5 times larger than the ones of the supernovae. Moreover, many 

long GRBs occurs without the presence of a supernova. It is concluded that in no way a 
GRB can originate from a supernova. The precise theoretical understanding of the GRB 
luminosity we present evidence, in both these systems, the existence of an independent 
component in the X-ray emission, usually interpreted in the current literature as part of 
the GRB afterglow. This component has been observed by Chandra and XMM to have 
a strong decay on scale of months. We have named here these two sources respectively 
URCA-1 and URCA-2, in honor of the work that George Gamow and Mario Shoenberg 
did in 1939 in this town of Urea identifying the basic mechanism, the Urea processes, 
leading to the process of gravitational collapse and the formation of a neutron star and 
a supernova. The further hypothesis is considered to relate this X-ray source to a neu¬ 
tron star, newly born in the Supernova. This hypothesis should be submitted to further 
theoretical and observational investigation. Some theoretical developments to clarify the 
astrophysical origin of this new scenario are outlined. 


1. Introduction 

In the last century the fundamental discoveries of nuclear physics have led to the 
understanding of the thermonuclear energy source of main sequence stars and ex¬ 
plained the basic physical processes underlying the solar luminosity (see e.g. M. 

Schwarzschild ( 3 )). 

The discovery of pulsars in 1968 (see Hewish et al. (.2)) led to the first evidence 
for the existence of neutron stars, first described in terms of theoretical physics by 
George Gamow as far back as 1936 0 ). It became clear that the pulsed luminosity 
of pulsars, at times 10 2 10 3 larger than solar luminosity, was not related to nuclear 

burning and could be simply explained in term of the loss of rotational energy of 
a neutron star (Gold Jj; y)). For the first time it became so clear the possible 
relevance of strong gravitational fields in the energetics of an astrophysical system. 

The birth of X-ray astronomy thanks to Riccardo Giacconi and his group (see 
e.g. Giacconi and Ruffini ( 0 )) led to a still different energy source, originating from 
the accretion of matter onto a star which has undergone a complete gravitational 
collapse process: a black hole (see e.g. Ruffini & Wheeler (?)). In this case, the 
energetics is dominated by the radiation emitted in the accretion process of mat¬ 
ter around an already formed black hole. Luminosities up to 10 4 times the solar 
luminosity, much larger then the ones of pulsars, could be explained by the release 
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of energy in matter accreting in the deep potential well of a black hole (Leach and 
Ruffini (J3)). This allowed to probe for the first time the structure of circular orbits 
around a black hole computed by Ruffini and Wheeler (see e.g. Landau and Lif- 
shitz |g)). This result was well illustrated by the theoretical interpretation of the 
observations of Cygnus-Xl, obtained by the Uhuru satellite and by the optical and 
radio telescopes on the ground (see Fig. [Q. 


The identification of the 
first black hole in our 
galaxy: Cygnus X-l 


• O = KF 7 erg/s = 10 4 /. 0 

= 0.01(</m/<f/) a(( c 2 

* Absence* of pulsation due to 
uniqueness of Kerr-Newmann 
black holes 

. M > 3.2 M 0 

Leach & Ruffini. 1973 



Figure 1. Cygnus X-l offered the possibility of identifying the first black hole in our galaxy (Leach 
and Ruffini J^)). The luminosity <f> of 10 4 solar luminosities points to the accretion process into a 
neutron star or a black hole as the energy source. The absence of pulsation is naturally explained 
either by a non-magnetized neutron star or a Kerr-Newmann black hole, which has necessarily to 
be axially symmetric. What identifies the black hole unambiguously is that the mass of Cygnus 
X-l, larger than 9Mq, exceeds the absolute upper limit of the neutron star mass, estimated at 
3.2M© by Rhoades and Ruffini ilOh . 


The discovery of gamma-ray bursts (GRBs) sign a further decisive progress. 
The GRBs give the first opportunity to probe and observe a yet different form of 
energy: the extractable energy of the black hole introduced in 1971 (Christodoulou 
and Ruffini (11)), which we shall refer in the following as the blackholic energy a . 


a This name is the English translation of the Italian words “energia buconerale”, introduced by 
Iacopo Ruffini, December 2004, here quoted by his kind permission. 








February 2, 2008 17:6 Proceedings Trim Size: 9.75in x 6.5in 


Rufiini'mglO 


4 

The blackholic energy, expected to be emitted during the dynamical process of 
gravitational collapse leading to the formation of the black hole, generates X-ray 
luminosities 10 21 times larger than the solar luminosity, although lasting for a very 
short time. 

The extreme regimes of GRBs evidence new and unexplored regimes of theoreti¬ 
cal physics. It is the aim of this talk to outline the progress achieved in understand¬ 
ing these astrophysical systems and the theoretically predicted regimes for the first 
time submitted to direct observational verification. 

It is a pleasure to present these results in this village of Urea, in the beautiful city 
of Rio. While sitting at the Casino de Urea, George Gamow and Mario Schoenberg 
in 1939 identified the basic process leading to the formation and cooling of a newly 
born neutron star (see |Appendix A| ) . They called this process essentially related 
to the emission of neutrinos and antineutrinos the Urea process. It is a welcomed 
coincidence that, in the last hours, while preparing this talk, examining the data of 
the recently observed GRB 030329, we have received a confirmation of a scenario 
we have recently outlined in three papers giving t he theoretical paradigms for the 
understanding of GRBs (Ruffini et al. {13 :llA Iblu. 

We have clear evidence, first advanced in the system GRB980425/SN1998bw 
(Ruffini et al. HH, Fraschetti et al. flldU and now confirmed also in the system 
GRB030329/SN1003dh, that there are in these systems three different components: 
1) the GRB source, generated by the collapse to a black hole, 2) the supernova, gen¬ 
erated by the collapse of an evolved star, 3) an additional X-ray source which is not 
related, unlike what is at times stated in the literature, to the GRB afterglow. In 
honour of the work done in the town of Urea by George Gamow and Mario Schoen¬ 
berg, identifying in the neutrino emission of the Urea process the basic mechanism 
leading to the process of gravitational collapse and the formation of a relativis¬ 
tic compact star, we named these two X-ray sources URCA-1, the one formed in 
the system GRB980425/SN1998bw, and URCA-2, the one formed in the system 
GRB030329/SN2003dh. We shall now recall some of the main steps in reaching 
this understanding out of the GRB phenomenon and explore possible explanation 
of the origin of these two sources. 


2. The energetics of gamma-ray bursts 

It is well known how GRBs were detected and studied for the first time using the 
Vela satellites, developed for military research to monitor the non-violation of the 
Limited Test Ban Treaty signed in 1963 (see e.g. Strong (1Z)). It was clear from 
the early data of these satellites, which were put at 150,000 miles from the surface 
of Earth, that the GRBs did not originate either on the Earth nor in the Solar 
System. 

The mystery of these sources became more profound as the observations of the 
BATSE instrument on board of the Compton Gamma-Ray Observatory (CGRO) 
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2704 BATSE Gamma-Ray Bursts 


+90 



Fluence, 50-300 keV (ergs cm' 2 ) 


Figure 2. Position in the sky, in galactic coordinates, of 2000 GRB events seen by the CGRO 
satellite. Their isotropy is evident. Reproduced from BATSE web site by their courtesy. 


satellite 13 over 9 years proved the isotropy of these sources in the sky (See Fig. 0. 
In addition to these data, the CGRO satellite gave an unprecedented number of 
details on the GRB structure, on their spectral properties and time variabilities 
which became encoded in the fourth BATSE catalog JTsIl (see e.g. Fig. 0. Out 
of the analysis of these BATSE sources it soon became clear (see e.g. Kouveliotou 
et al. Q, Tavani th e existence of two distinct families of sources: the long 
bursts, lasting more then one second and softer in spectra, and the short bursts (see 
Fig. 0 , harder in spectra (see Fig. 0 . We shall return shortly on this topic. 

The situation drastically changed with the discovery of the afterglow by the 
Italian-Dutch satellite BeppoSAX (Costa et al. (.21.)) and the possibility which led 
to the optical identification of the GRBs by the largest telescopes in the world, 
including the Hubble Space Telescope, the Keck Telescope in Hawaii and the VLT 
in Chile, and allowed as well the identification in the radio band of these sources. 
The outcome of this collaboration between complementary observational technique 
has led to the possibility of identifying in 1997 the distance of these sources from 
the Earth and their tremendous energy of the order up to 10 54 erg/second during 


b see http://cossc.gsfc.nasa.gov/batse/ 
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Trigger I 43 


Trigger 2067 


Trigger 2571 






Trigger 14D6 




Trigger 2514 



Figure 3. Some GRB light curves observed by the BATSE instrument on board of the CGRO 
satellite. 


the burst. It is interesting, as we will show in the following, that an energetics of 
this magnitude for the GRBs had previously been predicted out of first principles 
already in 1974 by Damour and Ruffini 171 ). 

The resonance between the X- and gamma ray astronomy from the satellites 
and the optical and radio astronomy from the ground, had already marked in the 
seventies the great success and development of the astrophysics of binary X-ray 
sources (see e.g. Giacconi & Ruffini p)). This resonance is re-proposed here for 
GRBs on a much larger scale. The use of much larger satellites, like Chandra and 
XMM-Newton, and dedicated space missions, like HETE-2 and, in the near future, 
Swift, and the very fortunate circumstance of the coming of age of the development 
of unprecedented optical technologies for the telescopes offers opportunities without 
precedence in the history of mankind. In parallel, the enormous scientific interest 
on the nature of GRB sources and the exploration, not only of new regimes, but 
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Figure 4. The energy fluence-averaged hardness ratio for short (T < 1 s)and long (T > 1 s) 
GRBs are represented. Reproduced, by his kind permission, from Tavani 1201) where the details 
are given. 


also of totally novel conceptual physical process of the blackholic energy, make the 
knowledge of GRBs an authentic new frontier in the scientific knowledge. 

3. The complexity an self-consistency of GRB modeling 

The study of GRBs is very likely “the” most extensive computational and theoretical 
investigation ever done in physics and astrophysics. There are at least three different 
fields of research which underlie the foundation of the theoretical understanding of 
GRBs. All three, for different reasons, are very difficult. 

The first field of research is the field of special relativity. As I always mention 
to my students in the course of theoretical physics, this held is paradoxically very 
difficult since it is extremely simple. In approaching special relativistic phenomena 
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Figure 5. On the upper right part of the figure are plotted the number of the observed GRBs 
as a function of their duration. The bimodal distribution corresponding respectively to the short 
bursts, upper left figure, and the long bursts, middle figure, is quite evident. The structure of the 
long bursts as part of the afterglow phenomena of GRBs is illustrated in section [8.11 


the extremely simple and clear procedures expressed by Einstein in his 1905 classic 
paper 0) are often ignored. Einstein makes use in his work of very few physical 
assumptions, an almost elementary mathematical framework and gives constant 
attention to a proper operational definition of all observable quantities. Those who 
work on GRBs use at times very intricate, complex and often wrong theoretical 
approaches lacking the necessary self-consistency. This is well demonstrated in the 
current literature on GRBs. 

The second field of research essential for understanding the energetics of GRBs 
deals with quantum electrodynamics and the relativistic process of pair creation in 
overcritical electromagnetic fields. This topic is also very difficult but for a quite dif¬ 
ferent conceptual reason: the process of pair creation, expressed in the classic works 
of Heisenberg-Euler-Schwinger 00 later developed by many others, is based on 
a very powerful theoretical framework but has never been verified by experimental 
data. The quest for creating electron-positron pairs by vacuum polarization pro¬ 
cesses in heavy ion collisions or in lasers has not yet been successfully achieved in 
Earth-bound experiments (see e.g. Ruffini, Vitagliano, Xue J 2 fil l'). As we will show 
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here, there is the tantalizing possibility of observing this phenomenon, for the first 
time, in the astrophysical setting of GRBs on a more grandiose scale. 

There is a third field which is essential for the understanding of the GRB phe¬ 
nomenon: general relativity. In this case, contrary to the case of special relativity, 
the field is indeed very difficult, since it is very difficult both from a conceptual, 
technical and mathematical point of view. The physical assumptions are indeed 
complex. The entire concept of geometrization of physics needs a new conceptual 
approach to the field. The mathematical complexity of the pseudo-Riemannian ge¬ 
ometry contrasts now with the simple structure of the pseudo-Euclidean Minkowski 
space. The operational definition of the observable quantities has to take into ac¬ 
count the intrinsic geometrical properties and also the cosmological settings of the 
source. With GRBs we have the possibility to follow, from a safe position in an 
asymptotically flat space at large distance, the formation of the horizon of a black 
hole with all the associated relativistic phenomena of light bending and time di¬ 
latation. Most important, as we will show in details in this presentation, general 
relativity in connection with quantum phenomena offers, with the blackholic energy, 
the explanation of the tremendous GRB energy sources. 

For these reasons GRBs offer an authentic new frontier in the field of physics and 
astrophysics. It is appropriate to mention some of the goals of such a new frontier 
in the above three fields. We recall in the special relativity field, for the first 
time, we observe phenomena occurring at Lorentz gamma factors of approximately 
300. In the field of relativistic quantum electro-dynamics we see for the first time 
the interchange between classical fields and the created quantum matter-antimatter 
pairs. In the field of general relativity also for the first time we can test the blackholic 
energy which is the basic energetic physical variable underlying the entire GRB 
phenomenon. 

The most appealing aspect of this work is that, if indeed these three different 
fields are treated and approached with the necessary technical and scientific matu¬ 
rity, the model which results has a very large redundancy built-in. The approach 
requires an unprecedented level of self-consistency. Any departures from the correct 
theoretical treatment in this very complex system lead to exponential departures 
from the correct solution and from the correct fit of the observations. 

It is so that, as the model is being properly developed and verified, its solution 
will have existence and uniqueness. 


3.1. GRBs and special relativity 

The ongoing dialogue between our work and the one of the workers on GRBs, rests 
still on some elementary considerations presented by Einstein in his classic article 
of 1905 fiS). These considerations are quite general and even precede Einstein’s 
derivation, out of first principles, of the Lorentz transformations. We recall here Ein¬ 
stein’s words: “We might, of course, content ourselves with time values determined 
by an observer stationed together with the watch at the origin of the co-ordinates, 






February 2, 2008 17:6 Proceedings Trim Size: 9.75in x 6.5in 


Rufiini'mglO 


10 

and co-ordinating the corresponding positions of the hands with light signals, given 
out by every event to be timed, and reaching him through empty space. But this 
co-ordination has the disadvantage that it is not independent of the standpoint of 
the observer with the watch or clock, as we know from experience”. 



Fig ure 6. Relation between the arrival time t a and the laboratory time t. Details in Ruffini et 

ai. mil. 


The message by Einstein is simply illustrated in Fig. El If we consider in an 
inertial frame a source (solid line) moving with high speed and emitting light signals 
(dashed lines) along the direction of its motion, a far away observer will measure a 
delay A t a between the arrival time of two signals emitted at the origin and after a 
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time interval At in the laboratory frame. The real velocity of the source is given 
by: 


V = 


and the apparent velocity is given by: 


u app 


A r 
At 


A r 

A T' 


(1) 

( 2 ) 


As pointed out by Einstein the adoption of coordinating light signals simply by 
their arrival time as in Eq. without an adequate definition of synchronization, 
is incorrect and leads to unsurmountable difficulties as well as to apparently “su¬ 
perluminal” velocities as soon as motions close to the speed of light are considered. 

The use of At 0 as a time coordinate, often tacitly adopted by astronomers, 
should be done, if at all, with proper care. The relation between A t a and the 
correct time parameterization in the laboratory frame has to be taken into account: 

A r 1 r t »+ At 

At a = At - = At - - v (t') dt'. (3) 

c c Jt 0 


In other words, the relation between the arrival time and the laboratory time cannot 
be done without a knowledge of the speed along the entire world-line of the source. 
In the case of GRBs, such a worldline starts at the moment of gravitational collapse. 
It is of course clear that the parameterization in the laboratory frame has to take 
into account the cosmological redshift z of the source. We then have, at the detector: 


A t d a = (!+*) At a 


( 4 ) 


In the current GRB literature, Eq. © has been systematically neglected by ad¬ 
dressing only the afterglow description neglecting the previous history of the source. 
Often the integral equation has been approximated by a clearly incorrect instanta¬ 
neous value: 


A t a 


At 
272 • 


( 5 ) 


The attitude has been adopted that it should be possible to consider separately 
the afterglow part of the GRB phenomenon, without the knowledge of the entire 
equation of motion of the source. 

This point of view has reached its most extreme expression in the works reviewed 
by Piran ©0), where the so-called “prompt radiation”, lasting on the order of 
10 2 s, is considered as a burst emitted by the prolonged activity of an “inner engine”. 
In these models, generally referred to as the “internal shock model”, the emission of 
the afterglow is assumed to follow the “prompt radiation” phase ( hit feEsEi- 

As we outline in the following, such an extreme point of view originates from 
the inability of obtaining the time scale of the “prompt radiation” from a burst 
structure. These authors consequently appeal to the existence of an “ad hoc” inner 
engine in the GRB source to solve this problem. 
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We show in the following how this difficulty has been overcome in our approach 
by interpreting the “prompt radiation” as an integral part of the afterglow and not 
as a burst. This explanation can be reached only through a relativistically correct 
theoretical description of the entire afterglow (see section^]). Within the framework 
of special relativity we show that it is not possible to describe a GRB phenomenon 
by disregarding the knowledge of the entire past worldline of the source. We show 
that at 10 2 seconds the emission occurs from a region of dimensions of approxi¬ 
mately 10 16 cm, well within the region of activity of the afterglow. This point was 
not appreciated in the current literature due to the neglect of the apparent superlu¬ 
minal effects implied by the use of the “pathological” parametrization of the GRB 
phenomenon by the arrival time of light signals. 

An additional difference between our treatment and the ones in the current lit¬ 
erature relates to the assumption of the existence of scaling laws in the afterglow 
phase: the power law dependence of the Lorentz gamma factor on the radial co¬ 
ordinate is usually systematically assumed. From the proper use of the relativistic 
transformations and by the direct numerical and analytic integration of the special 
relativistic equations of motion we demonstrate (see section 15.211 that no simple 
power-law relation can be derived for the equations of motion of the system. This 
situation is not new for workers in relativistic theories: scaling laws exist in the 
extreme ultrarelativistic regimes and in the Newtonian ones but not in the inter¬ 
mediate fully relativistic regimes (see e.g. Ruffini 0 ). 


3.2. GRBs and general relativity 

Three of the most important works in the field of general relativity have certainly 
been the discovery of the Kerr solution lifl) , its generalization to the charged case 
(Newman et al. (37)) and the formulation by Brandon Carter (3£| of the Hamilton- 
Jacobi equations for a charged test particle in the metric and electromagnetic field 
of a Kerr-Newman solution (see e.g. Landau and Lifshitz |90). The equations of 
motion, which are generally second order differential equations, were reduced by 
Carter to a set of first order differential equations which were then integrated by 
using an effective potential technique by Ruffini and Wheeler for the Kerr metric (see 
e.g. Landau and Lifshitz Q)) and by Ruffini for the Reissner-Nordstrom geometry 
(Ruffini (l;'i. r J) . see Fig. 0 . 

All the above mathematical results were essential for understanding the new 
physics of gravitationally collapsed objects and allowed the publication of a very 
popular article: “Introducing the black hole” (Ruffini and Wheeler Q)). In that 
paper, we advanced the ansatz that the most general black hole is a solution of 
the Einstein-Maxwell equations, asymptotically flat and with a regular horizon: the 
Kerr-Newman solution, characterized only by three parameters: the mass M, the 
charge Q and the angular momentum L. This ansatz of the “black hole uniqueness 
theorem” still today after thirty years presents challenges to the mathematical as¬ 
pects of its complete proof (see e.g. Carter © and Bini et al. (Sib). In addition 
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Figure 7. The effective potential corresponding to the circular orbits in the equatorial plane of a 
black hole is given as a function of the angular momentum of the test particle. This digram was 
originally derived by Ruffini and Wheeler (right picture). For details see Landau and Lifshitz 0 ) 
and Rees, Ruffini and Wheeler 


to these mathematical difficulties, in the field of physics this ansatz contains the 
most profound consequences. The fact that, among all the possible highly nonlinear 
terms characterizing the gravitationally collapsed objects, only the ones correspond¬ 
ing solely to the Einstein Maxwell equations survive the formation of the horizon 
has, indeed, extremely profound physical implications. Any departure from such 
a minimal configuration either collapses on the horizon or is radiated away during 
the collapse process. This ansatz is crucial in identifying precisely the process of 
gravitational collapse leading to the formation of the black hole and the emission 
of GRBs. Indeed, in this specific case, the Born-like nonlinear 62 ) term of the 
Heisenberg-Euler-Schwinger (24; 25) Lagrangian are radiated away prior to the for¬ 
mation of the horizon of the black hole (see e.g. Ruffini et al. f20|)). Only the 
nonlinearity corresponding solely to the classical Einstein-Maxwell theory is left as 
the outcome of the gravitational collapse process. 

The same effective potential technique (see Landau and Lifshitz ( 3 )), which 
allowed the analysis of circular orbits around the black hole, was crucial in reaching 
the equally interesting discovery of the reversible and irreversible transformations of 
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black holes by Christodoulou and Ruffini UJD. which in turn led to the mass-energy 
formula of the black hole: 


with 


where 



( 6 ) 

(7) 

( 8 ) 


is the horizon surface area, M j r is the irreducible mass, r+ is the horizon radius 
and p+ is the quasi-spheroidal cylindrical coordinate of the horizon evaluated at 
the equatorial plane. Extreme black holes satisfy the equality in Eq. 0 . 

From Eq. 0 follows that the total energy of the black hole Ebh can be split 
into three different parts: rest mass, Coulomb energy and rotational energy. In 
principle both Coulomb energy and rotational energy can be extracted from the 
black hole (Christodoulou and Ruffini 0 ) . The maximum extractable rotational 
energy is 29% and the maximum extractable Coulomb energy is 50% of the total 
energy, as clearly follows from the upper limit for the existence of a black hole, 
given by Eq. 0 . We refer in the following to both these extractable energies as the 
blackholic energy. 

The existence of the black hole and the basic correctness of the circular orbits 
has been proven by the observations of Cygnus-Xl (see e.g. Giacconi and Ruffini 
(6)). However, in binary X-ray sources, the black hole uniquely acts passively 
by generating the deep potential well in which the accretion process occurs. It 
has become tantalizing to look for astrophysical objects in order to verify the other 
fundamental prediction of general relativity that the blackholic energy is the largest 
energy extractable from any physical object. 

As we shall see in the next section, the feasibility of the extraction of the black¬ 
holic energy has been made possible by the quantum processes of creating, out of 
classical fields, a plasma of electron-positron pairs in the field of black holes. The 
manifestation of such process of energy extraction from the black hole is astrophys- 
ically manifested by the occurrence of GRBs. 


3.3. GRBs and quantum electro-dynamics 

That a static electromagnetic field stronger than a critical value: 


E, = 


m^c 3 

he 


(9) 


can polarize the vacuum and create electron-positron pairs was clearly evidenced 
by Heisenberg and Euler ji^). The major effort in verifying the correctness of this 
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theoretical prediction has been directed in the analysis of heavy ion collisions (see 
Ruffini et al. 0 ) and references therein). From an order-of-magnitude estimate, it 
appears that around a nucleus with a charge: 

hr 

137 (10) 

e z 

the electric field can be stronger than the electric field polarizing the vacuum. A 
more accurate detailed analysis taking into account the bound states levels around 
a nucleus brings to a value of 

Z c ~ 173 (11) 


for the nuclear charge leading to the existence of a critical field. From the Heisenberg 
uncertainty principle it follows that, in order to create a pair, the existence of the 
critical field should last a time 

h 


At ■ 


~ 10 


-18 , 


( 12 ) 


m e c * 

which is much longer then the typical confinement time in heavy ion collisions which 
is 

h 


At ■ 


rripC 


~ 10 -21 s. 


(13) 


This is certainly a reason why no evidence for pair creation in heavy ion collisions 
has been obtained although remarkable effort has been spent in various acceler¬ 
ators worldwide. Similar experiments involving laser beams encounter analogous 
difficulties (see e.g. Ruffini et al. ([26) and references therein). 

The alternative idea was advanced in 1975 |2^) that the critical field condition 
given in Eq. © could be reached easily, and for a time much larger than the one given 
by Eq. o, in the field of a Kerr-Newman black hole in a range of masses 3.2 Mq < 
Mbh < 7.2 x 1O 6 M 0 . In that paper we have generalized to the curved Kerr-Newman 
geometry the fundamental theoretical framework developed in Minkowski space by 
Heisenberg-Euler :{m ) and Schwinger (2£). This result was made possible by the 
work on the structure of the Kerr-Newman spacetime previously done by Carter 
l.isli and by the remarkable mathematical craftsmanship of Thibault Damour then 
working with me as a post-doc in Princeton. 

The maximum energy extractable in such a process of creating a vast amount 
of electron-positron pairs around a black hole is given by: 

E max = 1.8 x 10 54 ( Mbh/M q ) erg. (14) 


We concluded in that paper that such a process “naturally leads to a most simple 
model for the explanation of the recently discovered y-rays bursts”. 

At that time, GRBs had not yet been optically identified and nothing was known 
about their distance and consequently about their energetics. Literally thousands of 
theories existed in order to explain them and it was impossible to establish a rational 
dialogue with such an enormous number of alternative theories. We did not pursue 
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further our model until the results of the BeppoSAX mission, which clearly pointed 
to the cosmological origin of GRBs, implying for the typical magnitude of their 
energy precisely the one predicted by our model. 

It is interesting that the idea of using an electron-positron plasma as a basis 
of a GRB model was independently introduced years later in a set of papers by 
Cavallo and Rees (Thl - ! . Cavallo and Horstman 13 and Horstman and Cavallo @. 
These authors did not address the issue of the physical origin of their energy source. 
They reach their conclusions considering the pair creation and annihilation process 
occurring in the confinement of a large amount of energy in a region of dimension 
~ 10 km typical of a neutron star. No relation to the physics of black holes nor to 
the energy extraction process from a black hole was envisaged in their interesting 
considerations, mainly directed to the study of the opacity and the consequent 
dynamics of such an electron-positron plasma. 


A r = r - r 



e + e~ plasma 


Figure 8. The dyadosphere is comprised between the horizon radius and the radius of the dya- 
dosphere. All this reg ion is filled with electron-positron pairs and photons in thermal equilibrium. 
Details in Ruffini (l46i) . Preparata et al. l47jh Ruffini et al. J48fi . 

After the discovery of the afterglows and the optical identification of GRBs 
at cosmological distances, implying exactly the energetics predicted in Eq. m, 
we returned to the analysis of the vacuum polarization process around a black 
hole and precisely identified the region around the black hole in which the vacuum 
polarization process and the consequent creation of electron-positron pairs occur. 
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We defined this region, using the Greek name dyad for pairs (dvas, SvaScx 7 ), to be 
the “dyadosphere” of the black hole, bounded by the black hole horizon and the 
dyadosphere radius r ds given by (see Ruffini 0 , Preparata et al. 0 and Fig |Hl : 


rds = 



1.12- 10 8 y/Ji£cm, (15) 


where we have introduced the dimensionless mass and charge parameters fj, = 
Mbh/Mq, £ = Q/{MbhVG ) < 1. 

The analysis of the dyadosphere was developed, at that time, around an already 
formed black hole. In recent months we have been developing the dynamical forma¬ 
tion of the black hole and correspondingly of the dyadosphere during the process 
of gravitational collapse, reaching some specific signatures which may be detectable 
in the structure of the short and long GRBs (Cherubini et al. (49), Ruffini and 
Vitagliano Jjtllfilll. Ruffini et al. (0j;ISE3))- 


4. The dynamical phases following the dyadosphere formation 

Many details of this topic have been presented in great details in Ruffini et al. ( 27). 

After the vacuum polarization process around a black hole, one of the topics 
of the greatest scientific interest is the analysis of the dynamics of the electron- 
positron plasma formed in the dyadosphere. This issue was addressed by us in a very 
effective collaboration with Jim Wilson at Livermore. The numerical simulations 
of this problem were developed at Livermore, while the semi-analytic approach was 
developed in Rome (Ruffini et al. 0 ). 

The corresponding treatment in the framework of the Cavallo et al. analysis 
was performed by Piran et al. 0 also using a numerical approach, by Bisnovaty- 
Kogan and Murzina 0) using an analytic approach and by Meszaros, Laguna and 
Rees (.(3) using a numerical and semi-analytic approach. 

Although some analogies exists between these treatments, they are significantly 
different in the theoretical details and in the final results. Since the final result of 
the GRB model is extremely sensitive to any departure from the correct treatment, 
it is indeed very important to detect at every step the appearance of possible fatal 
errors. 

A conclusion common to all these treatments is that the electron-positron plasma 
is initially optically thick and expands till transparency reaching very high values 
of the Lorentz gamma factor. A second point, which is common, is the discovery of 
a new clear feature: the plasma shell expands but the Lorentz contraction is such 
that its width in the laboratory frame appears to be constant. 

There is however a major difference between our approach and the ones of Piran, 
Meszaros and Rees, in that the dyadosphere is assumed by us to be filled uniquely 
with an electron-positron plasma. Such a plasma expands in substantial agreement 
with the results presented in the work of Bisnovati-Kogan and Murzina 0 . In 
our model the pulse of electron-positron pairs and photons (PEM Pulse, see Ruffini 
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et al. ( 0 )) evolves and at a radius on the order of 10 10 cm it encounters the 
remnant of the star progenitor of the newly formed black hole. The PEM pulse 
is then loaded with baryons. A new pulse is formed of electron-positron-photons 
and baryons (PEMB Pulse, see Rufhni et al. (fjjf )) which expands all the way until 
transparency is reached. At transparency the emitted photons give origin to what 
we define as the Proper-GRB (see Rufhni et al. 0 ) and Fig. n. 


p r.RR 
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Figure 9. The optically thick phase of our model are qualitatively represented in this diagram. 
There are clearly recognizable 1) the PEM pulse phase, 2) the impact on the baryonic remnant, 
3) the PEMB pulse phase and the final app roach to transparency with the emission of the P-GRB 
(see Fig. 1101 . Details in Ruffini et al. l2ib . 


In our approach, the baryon loading is measured by a dimensionless quantity 

M b c 2 


B = 


E, 


dya 


(16) 


which gives direct information about the mass Mb of the remnant. The correspond¬ 
ing treatment done by Piran and collaborators (Sherni & Piran (59), Piran et al. 
0) and by Meszaros, Laguna and Rees 0 differs in one important respect: the 
baryonic loading is assumed to occur since the beginning of the electron-positron 
pair formation and no relation to the mass of the remnant of the collapsed progenitor 
star is attributed to it. 

A marked difference also exists between our description of the rate equation 
for the electron-positron pairs and the ones by those authors. While our results 
are comparable with the ones obtained by Piran under the same initial conditions, 
the set of approximations adopted by Meszaros, Laguna and Rees 0 appears to 
be too radical and leads to very different results violating energy and momentum 
conservation (see Bianco et al. (§Q)). 

From our analysis (Ruffini et al. |&§f)) it also becomes clear that such expanding 
dynamical evolution can only occur for values of B < 10 -2 . This prediction, as we 
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will show shortly in the three GRB sources considered here, is very satisfactorily 
confirmed by observations. 

From the value of the B parameter, related to the mass of the remnant, it 
therefore follows that the collapse to a black hole leading to a GRB is drastically 
different from the collapse to a neutron star. While in the case of a neutron star 
collapse a very large amount of matter is expelled, in many instances well above 
the mass of the neutron star itself, in the case of black holes leading to a GRB 
only a very small fraction of the initial mass (~ 10 ~ 2 or less) is expelled. The 
collapse to a black hole giving rise to a GRB appears to be much smoother than 
any collapse process considered until today: almost 99.9% of the star has to be 
collapsing simultaneously! 

We summarize in Figs. 0 blithe optically thick phase of GRBs in our model: we 
start from a given dyadosphere of energy Ed ya \ the pair-electromagnetic pulse (PEM 
pulse) self-accelerates outward typically reaching Lorentz gamma factors 7 ~ 200 
at r ~ 10 10 cm; at this point the collision of the PEM pulse with the remnant 
of the progenitor star occurs with an abrupt decrease in the value of the Lorentz 
gamma factor; a new pair-electromagnetic-baryon pulse (PEMB pulse) is formed 
which self-accelerates outward until the system becomes transparent. 



Figure 10. The P-GRB emitted at the transparency point at a time of arri val which has been 
computed following the prescriptions of Eq. 0. Details in Ruffini et al. SEl 









February 2, 2008 17:6 Proceedings Trim Size: 9.75in x 6.5in 


Rufiini'mglO 


20 


The photon emission at this transparency point is the Proper-GRB (P-GRB). 
An accelerated beam of baryons with an initial Lorentz gamma factor y 0 starts to 
interact with the interstellar medium at typical distances from the black hole of 
r 0 ~ 10 14 cm and at a photon arrival time at the detector on the Earth surface of 
if ~ 0.1 s. These values determine the initial conditions of the afterglow. 

5. The description of the afterglow 

After reaching transparency and the emission of the P-GRB, the accelerated bary- 
onic matter (the ABM pulse) interacts with the interstellar medium (ISM) and gives 
rise to the afterglow (see Fig. CD- Also in the descriptions of this last phase many 
differences exist between our treatment and the other ones in the current literature. 


E-APC 



v 
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Figure 11. The GRB afterglow phase is here represented together with the optically thick phase 
(see Fig. |9j. The value of the Lorentz gamma factor is here given from the transparency point all 
the way to the ultrarelativisitc, relativistic and non relativistic regimes. Details in Ruffini et al. 

Izzb. 


5.1. The initial value problem 

The initial conditions of the afterglow era are determined at the end of the opti¬ 
cally thick era when the P-GRB is emitted. As recalled in the last section, the 
transparency condition is determined by a time of arrival if, a value of the gamma 
Lorentz factor 70 , a value of the radial coordinate r 0 , an amount of baryonic matter 
Mb which are only functions of the two parameters Ed ya and B (see En. H1 tijl l. It 
is appropriate here to emphasize again that, in order to have the expansion leading 
to an observed GRB, one must have B < 10 -2 . 

This connection to the optically thick era is missing in the current approach in 
the literature which attributes the origin of the “prompt radiation” to an unspecified 
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inner engine activity (see Piran (28) and references therein). The initial conditions 
at the beginning of the afterglow era are obtained by a best fit of the later parts of 
the afterglow. This approach is quite unsatisfactory since, as we will explicitly show, 
the theoretical treatments currently adopted in the description of the afterglow are 
not correct. The fit using an incorrect theoretical treatment leads necessarily to 
the wrong conclusions as well as, in turn, to the determination of incorrect initial 
conditions. 


5.2. The equations of the afterglow dynamics 

Let us first summarize the commonalities between our approach and the ones in the 
current literature. In both cases (see Piran (jig), Chiang & Dermer 0 and Ruffini 
et al. 0) a thin shell approximation is used to describe the collision between the 
ABM pulse and the ISM: 

dE int = (7 - 1) dM ism c 2 , (17a) 

d'y = - 2 ^dM ism , (17b) 

dM = dE int + dM ism , (17c) 

dMi sm = 4:nm p n ism r 2 dr , (17d) 


where E- m t, 7 and M are respectively the internal energy, the Lorentz factor and 
the mass-energy of the expanding pulse, ? 7 j sm is the ISM number density which is 
assumed to be constant, rn p is the proton mass, £ is the emitted fraction of the 
energy developed in the collision with the ISM and M; sm is the amount of ISM 
mass swept up within the radius r: M; sm = ( 4 / 3 ) 7 r (r 3 — ro 3 )m p ni sm , where r 0 is 
the starting radius of the shock front. In general, an additional condition is needed 
in order to determine £ as a function of the radial coordinate. In the following, £ 
is assumed to be constant and such an approximation appears to be correct in the 
GRB context. 

In both our work and in the current literature (see Piran li^ l. Chiang & Dermer 
0 and Ruffini et al. |13)) a first integral of these equations has been found, leading 
to expressions for the Lorentz gamma factor as a function of the radial coordinate. 
In the “fully adiabatic condition” (i.e. e = 0) we have: 


, 2 = I 2 + 2 7o (M ism /M B ) + (Mj sm /M B ) 2 
1 + 2y 0 (M\ sm /M B ) + (Mi Sm /M B ) 


(18) 


while in the “fully radiative condition” (i.e. £ = 1 ) we have: 

1 + (M ism /M B ) (1 + 7 - 1 ) [1 + (1/2) (M ism /M B )\ 

"7 = ----— - (19) 

7o _1 + (M ism /M B ) (1 + 70 - 1 ) [1 + (1/2) (M ism /M B )] 

where 70 and M B are respectively the values of the Lorentz gamma factor and of 
the mass of the accelerated baryons at the beginning of the afterglow phase and r 0 
is the value of the radius r at the beginning of the afterglow phase. 
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A major difference between our treatment and the other ones in the current 
literature is that we have integrated the above equations analytically, obtaining 
the explicit analytic form of the equations of motion for the expanding shell in 
the afterglow for a constant ISM density. For the fully radiative case we have 
explicitly integrated the differential equation for r (t) in Eq. Q Dll . recalling that 
7~ 2 = 1— [dr/ (cdt)} , where t is the time in the laboratory reference frame. We have 
then obtained a new explicit analytic solution of the equations of motion for the 
relativistic shell in the entire range from the ultra-relativistic to the non-relativistic 
regimes: 


t = 


Mb —m'. 


t Q 


l(r-r 1| roVe ]n f M+0Wo)] 3 ( A 3 + l) } 
1 + l2cm°A 3 1 y ^3 + (r/r 0 ) 3 J(T+l) 3 J 

i0 


8 cVC 


r Q V3C 

6cm? A 2 


, 2 (r/r 0 )—A , 2—A 

arCtan Av 4 ~ arCtan T73 


( 20 ) 


where A = {/(M B - m°) /m°, C = M B 2 ( 7o - l)/(7o + 1) and m° = 

(4/3) Trm p n- Ism rl. 

Correspondingly, in the adiabatic case we have: 


t= 7o 


m° \ r — 


Mb ) c^tFi 4Msr o cvtFi 


+ to 


( 21 ) 


In the current literature, following Blandford and McKee |Q2J), a so-called “ultra- 
relativistic” approximation y 0 7 1 has been widely adopted by many authors 

to solve Eqs. lfTTll (see e.g. Sari liiil |2J), Waxman Rees & Meszaros (Q§), 
Granot et al. (S3), Panaitescu & Meszaros ^3), Phan | 2 ^). Gruzinov & Waxman 
Jfiot i. van Paradijs et al. 0 , Meszaros 0 and references therein). This leads to 
simple constant-index power-law relations: 


7 oc r a , ( 22 a) 

with a = 3 in the fully radiative case and a = 3/2 in the fully adiabatic case. 
This simple relation is in stark contrast to the complexity of Eq. (pi and Eq. m 
respectively. In the same spirit, instead of Eq. pi and Eq. GU, some authors have 
assumed the following much simpler approximation for the relation between the 
time and the radial coordinate of the expanding shell, both in the fully radiative 
and in the fully adiabatic cases: 


ct = r , ( 22 b) 

while others, like e.g. Panaitescu & Meszaros {is), have integrated the approximate 
Ea. (l22aD . obtaining: 


ct = r 


1 + (4a + 2) 1 7 2 (r) 


( 22 c) 


Again, it is appropriate here to emphasize the stark contrast between 
Ens. (l22bl) . H22cll and the exact analytic solutions of Eqs. O- expressed in 
Eas.lt2()l.lt2Tl. 
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5.3. The equitemporal surfaces (EQTSs) 

As pointed out long ago by Couderc (Jr2;), in all relativistic expansion the crucial 
geometrical quantities with respect to a physical observer are the “equitemporal 
surfaces” (EQTSs), namely the locus of source points of the signals arriving at the 
observer at the same time. 

For a relativistically expanding spherically symmetric source the EQTSs are 
surfaces of revolution about the line of sight. The general expression for their 
profile, in the form d = d(r), corresponding to an arrival time t a of the photons at 
the detector, can be obtained from (see e.g. Ruffini et al. lElf . Bianco and RufRni 

El El and Figs. EMU): 

ct a = ct (r) — r cos d + r* , (23) 


where r* is the initial size of the expanding source, d is the angle between the radial 
expansion velocity of a point on its surface and the line of sight, and t = t{r) is its 
equation of motion, expressed in the laboratory frame, obtained by the integration 
of Eos. 11711 . From the definition of the Lorentz gamma factor y -2 = 1 — ( dr/cdt ) 2 , 
we have in fact: 

ct (r) = f [l — y -2 (r')] 1 dr 1 , (24) 

J 0 

where y(r) comes from the integration of Ens. TI 71) . 

We have obtained the expressions in the adiabatic case and in the fully radiative 
cases respectively (see Bianco and Ruffini |3)): 


cos d = 





cto 

r 


ct a 


r 


r 


7 ° - ( m°/M B ) 

1 



(25) 


cos d = 


M b - m° 

~w^ [T 


r a ) + 


m°r 0 
8 rVC 



- 1 


+ 

12 rm°A 2 


[A+{r/r 0 )f 

(A 3 + l) j 

1 cto ct a 

A 3 + (r/r 0 ) 3 

(A + l)\ 

r r 


(26) 


r * ToVSC 
r 6 rm°A 2 


arctan 


2 ( 77 V 0 ) - A 
AV3 


2- A 

— arctan- -= 

AV3 J 


The two EQTSs are represented at selected values of the arrival time t a in Fig. ITUl 
where the illustrative case of GRB 991216 has been used as a prototype. The 
initial conditions at the beginning of the afterglow era are in this case given by 
7 o = 310.131, r 0 = 1.943 x 10 14 cm, t a = 6.481 x 10 3 s, r* = 2.354 x 10 8 cm (see 

Ruffini et al. 13 El El El))- 
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Figure 12. Not all values of rl are allowed. Only photons emitted at an angle such that cos x) > 
(v/c) can be viewed by the observer. Thus the maximum allowed ii value Umax corresponds to 
cosi) maI = (v/c). In this figure we show Umax (i.e. the angular amplitude of the visible area of 
the ABM pulse) in degrees as a function of the arrival time at the detector for the photons emitted 
along the line of sight (see text). In the earliest GRB phases v ~ c and so it max ~ 0. On the other 
hand, in the latest phases of the afterglow the ABM pulse vel ocity decreases and Umax tends to 
the maximum possible value, i.e. 90°. Details in Ruffini et al. 1 7ol: 2T l 


5.4. The bolometric luminosity of the source 


We assume that the internal energy due to kinetic collision is instantly radiated 
away and that the corresponding emission is isotropic. As in section FT21 let Ae be 
the internal energy density developed in the collision. In the comoving frame the 
energy per unit of volume and per solid angle is simply 

dE 
dV dfl 



due to the fact that the emission is isotropic in this frame. The total number of 
photons emitted is an invariant quantity independent of the frame used. Thus we 
can compute this quantity as seen by an observer in the comoving frame (which we 
denote with the subscript “o” ) and by an observer in the laboratory frame (which 
we denote with no subscripts). Doing this we find: 


dNj 


dN~ 


dtdfldD 


dtdfldY, 


A 3 cos d , 


( 28 ) 













February 2, 2008 17:6 Proceedings Trim Size: 9.75in x 6.5in 


Ruffini’mglO 


25 



Figure 13. The diameter of the visible area is represented as a function of the ABM pulse radius. 
In the earliest expansion phases (7 ~ 310) flmax is very small (see left pane and Fig. 1141 . so the 
visible area is just a small fraction of the total ABM pulse surface. On the other hand, in the final 
expansion ph ases fim ax —► 90° and almost all the ABM pulse surface becomes visible. Details in 
Ruffini et al. Uli 


where cost? comes from the projection of the elementary surface of the shell on the 
direction of propagation and A = 7(1 — (3 cos d) is the Doppler factor introduced in 
the two following differential transformation 

<m o =d£lx A -2 (29) 

for the solid angle transformation and 

dt 0 = dt x A -1 (30) 

for the time transformation. The integration in is performed over the visible 
area of the ABM pulse at laboratory time t, namely with 0 < d < d max and d max 
defined in section 15.31 (see Figs. Elm. An extra A factor comes from the energy 
transformation: 


E 0 = E x A. 

See also Chiang and Dermer Q. Thus finally we obtain: 


dE 


dE 


(31) 


dtdfldT, 


dtdfldT, 


A 4 cos d . 


(32) 
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Figure 14. This figure shows the temporal evolution of the visible area of the ABM pulse. The 
dashed half-circles are the expanding ABM pulse at radii corresponding to different laboratory 
times. The black curve marks the boundary of the visible region. The black hole is located at 
position (0,0) in this plot. Again, in the earliest GRB phases the visible region is squeezed along 
the line of sight, while in the final part of the afterglow phase almost all the emitted photons reach 
the observer. This time evolution of the visible area is crucial to the explanation of the GRB 
temporal structure. Details in Ruffini et al. 


Doing this we clearly identify ( dtdndi: ) 0 as ener gy density in the comoving frame 
up to a factor -f- (see Eq. (E3). Then we have: 


dE 

dtdVl 



—— v cos ■& A 4 

47T 


dE, 


(33) 


where the integration in dE is performed over the ABM pulse visible area at labo¬ 
ratory time t, namely with 0 < d < dmax and d maa; defined in section 15.31 Ea. ll33l) 
gives us the energy emitted toward the observer per unit solid angle and per unit 
laboratory time t in the laboratory frame. 

What we really need is the energy emitted per unit solid angle and per unit 
detector arrival time if, so we must use the complete relation between if and t 
given in Eq. EH>. First we have to multiply the integrand in Eq. (|33|) by the factor 
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Figure 15. Due to the extremely high and extremely varying Lorentz gamma factor, photons 
reaching the detector on the Earth at the same arrival time are actually emitted at very different 
times and positions. We represent here the surfaces of photon emission corresponding to selected 
values of the photon arrival time at the detector: the equitemporal surfaces (EQTS). Such surfaces 
differ from the ellipsoids described by Rees in the context of the expanding radio sources with 
typical Lorentz factor 7 ~ 4 and constant. In fact, in GRB 991216 the Lorentz gamma factor 
ranges from 310 to 1. The EQTSs represented here (solid lines) correspond respectively to values 
of the arrival time ranging from 5 s (the smallest surface on the left of the plot) to 60 s (the largest 
one on the right). Each surface differs from the previous one by 5 s. To each EQTS contributes 
emission processes occurring at different values of the Lorentz gamma factor. The dashed lines are 
the boundaries of the visible area of the ABM pulse and the black hole is located at position (0, 0) 
in this plot. Note the different scales on the two axes, indicating the very high EQTS “effective 
eccentricity”. The time interval from 5 s to 60 s has been chosen to enco mpass the E-APE emission, 
ranging from 7 = 308.8 to 7 = 56.84. Details in Ruffini et al. SIH) 


( dt/dta ) to transform the energy density generated per unit of laboratory time 
t into the energy density generated per unit arrival time t Then we have to 
integrate with respect to over the equitemporal surface (EQTS, see section mil) 
of constant arrival time instead of the ABM pulse visible area at laboratory time 
t. The analog of Eq. (pi for the source luminosity in detector arrival time is then: 


dE 7 
dt^dfl 


I EQTS 


Ag a a -4 dt 

— v cosd A —r<i£ . 

4 tt dti 


( 34 ) 
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Figure 16. Comparison between EQTSs in the adiabatic regime (solid lines) and in the fully 
radiative regime (dashed lines). The upper plot shows the EQTSs for t a = 5 s, t a = 15 s, t a = 30 
s and t a = 45 s, respectively from the inner to the outer one. The lower plot shows the EQTS at 
an arrival time of 2 days. Details in Bianco and Ruffini El- 
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It is important to note that, in the present case of GRB 991216, the Doppler factor 
A~ 4 in Eq. m enhances the apparent luminosity of the burst, as compared to the 
intrinsic luminosity, by a factor which at the peak of the afterglow is in the range 
between 10 10 and 10 12 ! 



Figure 17. Bolometric luminosity of P-GRB and afterglow as a function o f th e arrival time. 
Details in Ruffini et al. El- Reproduced and adapted from Ruffini et al. with the kind 

permission of the publisher. 


We are now able to reproduce in Fig. II71 the general behavior of the luminosity 
starting from the P-GRB to the latest phases of the afterglow as a function of 
the arrival time. It is generally agreed that the GRB afterglow originates from an 
ultrarelativistic shell of baryons with an initial Lorentz factor y 0 ~ 200-300 
respect to the interstellar medium (see e.g. Ruffini et al. (2j|, Bianco & Ruffini 
and references therein). Using GRB 991216 as a prototype, in Ruffini et al. (0 
we have shown how from the time varying bolometric intensity of the afterglow it 
is possible to infer the average density ( rii sm ) = 1 particle/cm 3 of the InterStellar 
Medium (ISM) in a region of approximately 10 1 ' cm surrounding the black hole 
giving rise to the GRB phenomenon. 

It was shown in Ruffini et al. 0 that the theoretical interpretation of the 
intensity variations in the prompt phase in the afterglow implies the presence in 


with 



















February 2, 2008 17:6 Proceedings Trim Size: 9.75in x 6.5in 


Rufiini'mglO 


30 

the ISM of inhomogeneities of typical scale 10 15 cm. Such inhomogeneities were 
there represented for simplicity as spherically symmetric over-dense regions with 
{ n °ifm) — 10 2 ( n ism) separated by under-dense regions with (n^ m ) ~ 10 -2 (rii sm ) 
also of typical scale ~ 10 15 cm in order to keep ( rii srn ) constant. 

The summary of these general results are shown in Fig. I1SI where the P-GRB, 
the emission at the peak of the afterglow in relation to the “prompt emission” and 
the latest part of the afterglow are clearly identified for the source GRB 991216. 
Details in Ruffini et al. (27:). 


GRB991216 



Ruffini, Bianco. Chardonnct, Fraschctti, Xuc. ApJ, 555. LI 13,(2001) 
Ruffini. Bianco. Chardonnet, Fraschetti. Xue. ApJ, 581. L19. (2002) 
Ruffini, Bianco, Chardonnct, Fraschctti, Vilugliano, Xuc,‘‘Cosmology and Gra\ilalion**, AIP, (2003) 


Figure 18. The detailed features of GRB 991216 evidenced by our theoretical models are here 
reproduced. The P-GRB, the “prompt radiation” and what is generally called the afterglow. It 
is clear that the prompt radiation coincides with the extended afterglow peak emission (E-APE) 
and has been considered as a burst only as a consequence of the high noise threshold in the 
observations. Details in Ruffini et al. 17a 1271) . 


6. The theory of the luminosity in fixed energy bands and spectra 
of the afterglow 

Having obtained a general agreement between the observed luminosity variability 
and our treatment of the bolometric luminosity, we have further developed the 
model in order to explain 
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a) the details of the observed luminosity in fixed energy bands, which are the ones 
actually measured by the detectors on the satellites 

b) the instantaneous as well as the average spectral distribution in the entire after¬ 
glow and 

c) the observed hard to soft drift observed in GRB spectra. 

In order to do so we have developed (Ruffini et al. (j./EI) a more detailed theory 
of the structure of the shock front giving rise to the afterglow. We have modeled 
the interaction between the ultrarelativistic shell of baryons and the ISM by a shock 
front with three well-defined layers (see e.g. secs. 85-89, 135 of Landau & Lifshitz 
(7sjb ch. 2 and sec. 13-15 of Zel’dovich & Rayzer { 79 I) and sec. IV, 11 13 of Sedov 
l 8y)). From the back end to the leading edge of this shock front there is: 

a) A compressed high-temperature layer, of thickness A', in front of the relativistic 
baryonic shell, created by the accumulated material swept up in the ISM. 

b) A thin shock front, with a jump AT in the temperature which has been tradition¬ 
ally estimated in the comoving frame by the Rankine-Hugoniot adiabatic equations: 

AT ~ (3 / 16)m p 5v 2 / k ~ 1.5 x 10 11 [6v/(10 5 kms ~ 1 )] 2 K , (35) 

where Sv is the velocity jump, m v is the proton mass and k is Boltzmann’s constant. 
Of course such a treatment, valid for 7 ~ 1, has to be modified (see below) in our 
novel treatment for the 7 ~ 200 case relevant to GRBs. 

c) A pre-shock layer of ISM swept-up matter at much lower density and tempera¬ 
ture, both of which change abruptly at the thin shock front behind it. 

At larger distances ahead of the expanding fireball the ISM is at still smaller 
densities. The upper limit to the temperature jump at the thin shock front, given 
in Eq. ll35l) . is due to the transformation of kinetic energy to thermal energy, since 
the particle mean free path is assumed to be less then the thickness of the layer (a). 
The thermal emission of the observed X- and gamma ray radiation, which as seen 
from the observations reveals a high level of stability, is emitted in the above region 
(a) due to the sharp temperature gradient at the thin shock front described in the 
above region (b). 

The optical and radio emission comes in our model from the extended region 
(c). The description of such a region, unlike the sharp and well-defined tempera¬ 
ture gradient occurring in region (b), requires magnetohydrodynamic simulations 
of the evolution of the electron energy distribution of the synchrotron emission. 
Such analysis has been performed using 3-D Eulerian MHD codes for the particle 
acceleration models to produce the energy spectrum of cosmic rays at supernova en¬ 
velope fronts (see e.g. McKee and Cowie Isit h Tenorio-Tagle et al. Stone and 
Norman ®, Jun & Jones Isjl )!. Other challenges are the magnetic field and the 
instabilities. We mention two key phenomena: first, the importance of the develop¬ 
ment of Kelvin-Helmholtz and Rayleigh-Taylor instabilities ahead of the thin shock 
front. The second is the dual effect that the shock front has on the ISM initial mag¬ 
netic field, first through the compression of the swept-up matter containing the field 
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and secondly the amplification of the radial magnetic field component due to the 
Rayleigh-Taylor instability. Simulations of both effects (see e.g. Jun and Jones {id) 
and references therein), modeling the synchrotron radio emission for an expanding 
supernova shell at various initial magnetic field and ISM parameter values, shows 
for example that the presence of an initial tangential magnetic field component may 
essentially affect the resulting magnetic field configuration and hence the outgoing 
radio flux and spectrum. Among the additional effects to be taken into account are 
the initial inhomogeneity of the ISM and the contribution of magnetohydrodynamic 
turbulence. 

In our approach we focus uniquely on the X- and gamma ray radiation, which 
appears to be conceptually much simpler than the optical and radio emission. It 
is perfectly predictable by a set of constitutive equations (see next section), which 
leads to directly verifiable and very stable features in the spectral distribution of the 
observed GRB afterglows. In line with the observations of GRB 991216 and other 
GRB sources, we assume in the following that the X- and gamma ray luminosity 
represents approximately 90% of the energy flux of the afterglow, while the optical 
and radio emission represents only the remaining 10%. 

This approach differs significantly from the other ones in the current literature, 
where attempts are made to explain at once all the multi-wavelength emission in 
the radio, optical, X and gamma ray as coming from a common origin which is 
linked to boosted synchrotron emission. Such an approach has been shown to have 
a variety of difficulties (Ghirlanda et al. js3), Preece et al. lidi li and cannot anyway 
have the instantaneous variability needed to explain the structure in the “prompt 
radiation” in an external shock scenario, which is indeed confirmed by our model. 


6.1. The equations determining the luminosity in fixed energy 
bands 

Here the fundamental new assumption is adopted (see also Ruffini et al. 
the X- and gamma ray radiation during the entire afterglow phase has 
spectrum in the co-moving frame. The temperature is then given by: 

T s = [AE int / (4i r 2 Ar ff R)] 1/4 , (36) 

where AJ?i nt is the internal energy developed in the collision with the ISM in a time 
interval At in the co-moving frame, a is the Stefan-Boltzmann constant and 

n = A eff /A, (37) 

is the ratio between the “effective emitting area” of the afterglow and the surface 
area of radius r. In GRB 991216 such a factor is observed to be decreasing during 
the afterglow between: 3.01 x 10~ 8 >1Z> 5.01 x 10 -12 (Ruffini et al. (&Z)). 

The temperature in the comoving frame corresponding to the density distribu¬ 
tion described in Ruffini et al. ((75!) is shown in Fig. 1 1 91 


(.83)) that 
a thermal 
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Figure 19. The temperature in the comoving frame of the shock front corresponding to the 
density distribution with the six spikes A,B,C,D,E,F presented in Ruffini et al. 5 . The dashed line 
corresponds to an homogeneous distribution with rii srn = 1. Details in Ruffini et al. & 


We are now ready to evaluate the source luminosity in a given energy band. The 
source luminosity at a detector arrival time per unit solid angle dfl and in the 
energy band [ 17 , z/ 2 ] is given by (see Ruffini et al. ( 27; E?)): 


dE [ ^’ V2] 


A £ 


v cos d A 4 -jjiW (ui,f 2 ,T arr )dT. 


(38) 


where Ae = AE int /V is the energy density released in the interaction of the ABM 
pulse with the ISM inhomogeneities measured in the comoving frame, A = 7(1 — 
( v/c ) cosd) is the Doppler factor, W ( 17,17 ,T arr ) is an “effective weight” required to 
evaluate only the contributions in the energy band [ 17 , 17 ], dT is the surface element 
of the EQTS at detector arrival time on which the integration is performed (see 
also Ruffini et al. I ?- 1 ! )) and T arr is the observed temperature of the radiation 
emitted from dT,: 


Tarr = T./ [ 7 (1 - (v / c) COS d) (1 + Z )\ . (39) 

The “effective weight” W (z7,17 ,T arr ) is given by the ratio of the integral over 
the given energy band of a Planckian distribution at a temperature T arr to the total 
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integral aT^ rr : 



(40) 


where p(T arr ,v) is the Planckian distribution at temperature T c 


am 


p (T arr , V ) = ( 2/h 3 ) hv / (e^/( feT —) - l) 


(41) 


7. On the time integrated spectra and the hard-to-soft spectral 
transition 

We turn now to the much debated issue of the origin of the observed hard-to-soft 
spectral transition during the GRB observations (see e.g. Frontera et al. (88), 
Ghirlanda et al. |85), Piran (28), Piro et al. ( 0 )). We consider the instantaneous 
spectral distribution of the observed radiation for three different EQTSs: 

• ta = 10 s, in the early radiation phase near the peak of the luminosity, 

• = 1.45 x 10 5 s, in the last observation of the afterglow by the Chandra 
satellite, and 

• = 10 4 s, chosen in between the other two (see Fig. ED- 

The observed hard-to-soft spectral transition is then explained and traced back to: 

(1) a time decreasing temperature of the thermal spectrum measured in the 
comoving frame, 

(2) the GRB equations of motion, 

(3) the corresponding infinite set of relativistic transformations. 

A clear signature of our model is the existence of a common low-energy behavior of 
the instantaneous spectrum represented by a power-law with index a = +0.9. This 
prediction will be possibly verified in future observations. 

Starting from these instantaneous values, we integrate the spectra in arrival 
time obtaining what is usually fit in the literature by the “Band relation” (Band 
et al. l9(ik Indeed we find for our integrated spectra a low energy spectral index 
a = —1.05 and an high energy spectral index /3 < —16 when interpreted within 
the framework of a Band relation (see Fig. ED- This theoretical result can be 
submitted to a direct confrontation with the observations of GRB 991216 and, 
most importantly, the entire theoretical framework which we have developed can 
now be applied to any GRB source. The theoretical predictions on the luminosity 
in fixed energy bands so obtained can be then straightforwardly confronted with 
the observational data. 


8. The three paradigms for the interpretation of GRBs 

Having outlined the main features of our model and shown its application to GRB 
991216 used as a prototype, before addressing the two new sources which are going to 
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Figure 20. The instantaneous spectra of the radiation observed in GRB 991216 at three different 
EQTS respectively, from top to bottom, for t = 10 s, = 10 4 s and t ^ = 1.45 x 10 5 s. These 
diagrams have been computed assuming a constant (rii srn ) ~ 1 particle/cm 3 and c lear ly explains 
the often quoted hard-to-soft spectral evolution in GRBs. Details in Ruffini et al. E2). 


be the focus of this presentation, we recall the three paradigms for the interpretation 
of GRBs we had previously introduced. 

The first paradigm, the relative space-time transformation (RSTT) paradigm 
(Ruffini et al. 0) emphasizes the importance of a global analysis of the GRB 
phenomenon encompassing both the optically thick and the afterglow phases. Since 
all the data are received in the detector arrival time it is essential to know the 
equations of motion of all relativistic phases with 7 > 1 of the GRB sources in order 
to reconstruct the time coordinate in the laboratory frame, see Eq. 0 - Contrary 
to other phenomena in nonrelativistic physics or astrophysics, where every phase 
can be examined separately from the others, in the case of GRBs all the phases are 
inter-related by their signals received in arrival time t There is the need, in order 
to describe the physics of the source, to derive the laboratory time f as a function 
of the arrival time t „ along the entire past worldline of the source using Eq. 0 - 

The second paradigm, the interpretation of the burst structure (IBS) paradigm 
(Ruffini et al. (12|) covers three fundamental issues: 

a) the existence, in the general GRB, of two different components: the P-GRB and 
the afterglow related by precise equations determining their relative amplitude and 
temporal sequence (see Ruffini et al. (22)); 
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Figure 21. The time-integrated spectrum of the radiation observed in GRB 991216. The low 
energy part of the curve below 10 keV is fit by a power-law with index a = —1.05 and the high 
ener gy part above 500 keV is fit by a power-law with an index (3 < —16. Details in Ruffini et al. 


b) what in the literature has been addressed as the “prompt emission” and consid¬ 
ered as a burst, in our model is not a burst at all — instead it is just the emission 
from the peak of the afterglow (see Fig. CHI; 

c) the crucial role of the parameter B in determining the relative amplitude of the 
P-GRB to the afterglow and discriminating between the short and the long bursts 
(see Fig. El- Both short and long bursts arise from the same physical phenomena: 
the dyadosphere. The absence of baryonic matter in the remnant leads to the short 
bursts and no afterglow. The presence of baryonic matter with B < 10~ 2 leads 
to the afterglow and consequently to its peak emission which gives origin to the 
so-called long bursts. 

The third paradigm, the GRB-Supernova Time Sequence (GSTS) paradigm 
(Ruffini et al. 0 )) , deals with the relation of the GRB and the associated su¬ 
pernova process, and acquires a special meaning in relation to the sources GRB 
980425 and GRB 030329 as we will show in the following. 

We now shortly illustrate some consequences of these three paradigms. 
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Figure 22. The energy radiated in the P-GRB (the solid line) and in the afterglow (the dashed 
line), in units of the total energy of the dyadosphere ( Edya) 5 are plotted as functions of the B 
parameter. 


8.1. Long bursts are E-APEs 

The order of magnitude estimate usually quoted for the the characteristic time scale 
to be expected for a burst emitted by a GRB at the moment of transparency at the 
end of the expansion of the optically thick phase is given by r ~ GM/c 3 , which for 
a 10 Mg black hole will give ~ 10 -3 s. There are reasons today not to take seriously 
such an order of magnitude estimate (see e.g. Rufhni et al. |5J)). In any case this 
time is much shorter then the ones typically observed in “prompt radiation” of the 
long burstSjfrom a few seconds all the way to 10 2 s. In the current literature (see 
e.g. Piran (1281 ) and references therein), in order to explain the “prompt radiation” 
and overcome the above difficulty it has been generally assumed that its origin 
should be related to a prolonged “inner engine” activity preceding the afterglow 
which is not well identified. 

To us this explanation has always appeared logically inconsistent since there 
remain to be explained not one but two very different mechanisms, independent of 
each other, of similar and extremely large energetics. This approach has generated 
an additional very negative result: it has distracted everybody working in the field 
from the earlier very interesting work on the optically thick phase of GRBs. 

The way out of this dichotomy in our model is drastically different: 1) indeed the 
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optically thick phase exists, is crucial to the GRB phenomenon and terminates with 
a burst: the P-GRB; 2) the “prompt radiation” follows the P-GRB; 3) the “prompt 
radiation” is not a burst: it is actually the temporally extended peak emission of 
the afterglow (E-APE). The observed structures of the prompt radiation can all be 
traced back to inhomogeneities in the interstellar medium (see Fig. II81 and RufRni 
et al. ©). 


8.2. Short bursts are P-GRBs 

The fundamental diagram determining the relative intensity of the P-GRB and the 
afterglow as a function of the dimensionless parameter B has been shown in Fig. 
E3 The underlying machine generating the short and the long GRBs is identical: 
in both cases is the dyadosphere. The main difference relates to the amount of 
baryonic matter engulfed by the electron-positron plasma in their optically thick 
phase prior to transparency. In the limit of small B < 10 ~ 5 the intensity of the 
P-GRB is larger and dominates the afterglow. This corresponds to the short bursts. 
For 10 -5 < B < 10 -2 the afterglow dominates the GRBs and we have the so-called 
“long bursts”. For B > ICE 2 we may observe a third class of “bursts”, eventually 
related to a turbulent process occurring prior to transparency (Ruffini et al. 

This third family should be characterized by smaller values of the Lorentz gamma 
factors than in the case of the short or long bursts. 


8.3. The trigger of multiple gravitational collapses 

The relation between the GRBs and the supernovae is one of the most complex 
aspects to be addressed by our model, which needs the understanding of new fields 
of general relativistic physics in relation to yet unexplored many-body solutions in 
a substantially new astrophysical scenario. 

As we will show in the two systems GRB980425/SN1998bw 
and GRB030329/SN2003dh which we are going to discuss next, there is in each 
one the possibility of an astrophysical “triptych ” 0 formed by: 

1) the formation of the black hole and the emission of the GRB, 

2 ) the gravitational collapse of an evolved companion star, leading to a supernova, 

3) a clearly identified URCA source whose nature appears to be of the greatest 
interest. 

This new astrophysical scenario presents new challenges: 
a) The identification of the physical reasons of the instability leading to the gravita¬ 
tional collapse of a ~ IOMq star, giving origin to the black hole. Such an implosion 
must occur radially with negligible mass of the remnant (B < 10 “ 2 ). 


C A picture or carving in three panels side by side; esp: an altarpiece with a central panel and two 
flanking panels half its size that fold over it [Webster’s New collegiate dictionary, G. & C. Merriam 
Co. (Springfield, Massachussets, U.S.A., 1977)] 
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b) The identification of the physical reasons for the instability leading to the grav¬ 
itational collapse of an evolved companion star, giving origin to the supernova. 

c) The theoretical issues related to the URCA sources, which range today in many 
possible directions: from the physics of black holes, to the physical processes occur¬ 
ring in the expanding supernova remnants, and finally to the very exciting possibility 
that we are observing for the first time a newly born neutron star. The main effort 
in the next sections is to show that the detailed understanding we have reached for 
the GRB phenomenon and its afterglow allows us to state, convincingly, that the 
URCA source, contrary to what established in the current literature, is not part of 
the GRB nor of its afterglow. 

We will draw in the conclusions some considerations on the possible nature of 
the URCA sources. 

9. Applications 

We illustrate the application of our GRB model to two different systems, which are 
quite different in the energetics but are both related to supernovae: GRB 980425 
and GRB 030329. We will let the gradual theoretical understanding of the system 
to unveil the underlying astrophysical scenario. 


9.1. GRB 980425 / SN 1998bw 

Approaches in the current literature have always attempted to explain both the 
supernova and the GRB as two aspects of a single phenomenon assuming that the 
GRB takes his origin from a specially strong and yet unobserved supernova process: 
a hypernova (see Paczynski (f9lh . Kulkarni Il92lf . Iwamoto 0». 

We have taken a very different approach, following Cicero’s classic aphorism 
“divide et impera”, which was adopted as the motto of the Roman empire: “divide 
and conquer”. In this specific case of GRBs, which are indeed a very complex 
system, we plan to divide and identify the truly independent physical constituents 
and conquer the understanding of the underlying astrophysical process. As we will 
see, this approach will lead to an unexpected and much richer scenario. 

In addition to the source GRB 980425 and the supernova SN1998bw, two X- 
ray sources have been found by BeppoSAX in the error box for the location of 
GRB 980425: a source SI and a source S2 (Pian et al. libil b which have been 
traditionally interpreted either as a background source or as a part of the GRB 
afterglow. See Fig. 1231 Our approach has been: to first comprehend the entire 
afterglow of GRB 980425 within our theory. This allows the computation of the 
luminosity in given energy bands, the spectra, the Lorentz gamma factors, and 
more generally of all the dynamical aspects of the source. Having characterized 
the features of GRB 980425, we can gradually approach the remaining part of 
the scenario, disentangling the GRB observations from those of the supernova and 
then disentangling both the GRB and the supernova observations from those of the 
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sources SI and S2. This leads to a natural identification of distinct events and to 
their autonomous astrophysical characterization. 


GRB 980425 - SN1998bw 



Figure 23. The “divide et impera” concept applied to the system GRB 980425 / SN 1998bw. 
Four different components are identified: the GRB 980425, the SN 1998bw, and the two sources 
SI and S2. 


Our best fit for GRB 980425 corresponds to Ed ya = 1.1 x 10 48 ergs, B = 7 x 10“ 3 
and the ISM average density is found to be (ni sm ) = 0.02 particle/cm 3 . The plasma 
temperature and the total number of pairs in the dyadosphere are respectively 
T = 1.028 MeV and N e ± = 5.32 74 x 10 53 . The light curve of the GRB is shown in 
Figs. E3H23 The P-GRB is under the threshold and in the case of this source is 
not observable (see RufHni et al. J1 fll . Fraschetti et al. {Tf*)). 

The characteristic parameter 77, defining the filamentary structure of the ISM, 
monotonically decreases from 4.81 x 10 _1 ° to 2.65 x 10~ 12 ). The results are given 
in Fig. I20l where the bolometric luminosity is represented together with the optical 
data of SN1998bw, the source SI and the source S2. It is then clear that GRB 
980425 is separated both from the supernova data and from the sources SI and S2. 

While the occurrence of the supernova in relation to the GRB has already been 
discussed within the GRB-Supernova Time Sequence (GSTS) paradigm (Ruffini 
et al. 0 ) , we like to address here a different fundamental issue: the nature of 
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GRB 980425 - SN1998bw 



Figure 24. The bolometric luminosity as a function of the arrival time. The peak of the P-GRB 
is just below the observational noise level. 


the source SI which we have named, in celebration of the work of Gamow and 
Shoenberg, URCA-1. It is clear, from the theoretical predictions of the afterglow 
luminosity, that the URCA-1 cannot be part of the afterglow (see Figs. 124112611 . 
There are three different possibilities for the explanation of such source: 

1) Its possible relation to the black hole formed during the process of gravitational 
collapse leading to the GRB emission. 

2) Its possible relation to emission originating in the early phases of the expansion 
of the supernova remnant. 

3) The very exciting possibility that for the first time we are observing a newly born 
neutron star out of the supernova phenomenon. 

While some general considerations will be discussed in the conclusions, we would 
like to stress here the paramount importance of following the further time history of 
URCA-1 and of the source S2. If, as we propose, S2 is a background source, its flux 
should be practically constant in time and this source has nothing to do with the 
GRB 980425 / SN1998bw system. The drastic behavior of the URCA-1 luminosity 
reported in the talk by Elena Pian in this meeting, showing the latest URCA-1 
observations by the XMM and Chandra satellites, is crucial for the understanding 
of the nature of this source. Some very qualitative luminosity curves are sketched 
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The prompt emission of GRB 980425 
and the fit in the EMBH model 



Figure 25. The observation by BeppoSAX of the peak of the afterglow in the 40-700 keV energy 
band is fitted by our model. 


in Fig. no illustrating the possible time evolution of URCA-1. They are still very 
undetermined today due to a lack of attention to these observational data and, 
consequently, to the lack of a detailed theoretical model of the phenomenon. We 
therefore propose to have a dedicated attention to the astrophysical “triptych” GRB 
980425 / SN 1998bw / URCA-1. 

9.2. GRB 030329 / SN 2003dh 

We have adopted for our modeling of GRB030329 a spherically symmetric distribu¬ 
tion for the source and, as initial conditions at t = 10~ 21 s, an e + -e _ -photon neutral 
plasma lying between the radii rr = 2.9 x 10 6 cm and r 2 = 9.0 x 10 7 cm. The tem¬ 
perature of such a plasma is 2.1 MeV, the total energy E tot = 2.1 x 10 52 erg and the 
total number of pairs N e + e - = l.lxlO 57 . The baryonic matter component Mb is the 
second free parameter of our theory: B = 4.8 x 10 -3 . At the emission of the P-GRB, 
the Lorentz gamma factor is y 0 = 183.6 and the radial coordinate is r 0 = 5.3 x 10 13 
cm. The ISM average density is best fit by < n, sm >= 1 particle/cm 3 . The third 
free parameter of our theory is given by 1.1 x 10 ~ 7 < 1Z < 5.0 x 10 -11 . 

We then obtain (see also Bernardini et al. (££) in these proceedings) for the GRB 
030329 the luminosities in given energy bands, computed in the range 2-400 keV 
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GRB 980425 - SN1998bw - URCA-1 
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Figure 26. The bolometric light curves is reported as well as the BeppoSAX MECS observations 
in t he 2—10 keV band of SI and S2 (Pian et al. 1 94l ) ) and the optical data of SN 1998bw (Iwamoto 
$9j5i)). Here we also repo rt so me very qualitative curves to be expected for the URCA-1 luminosity. 
Details in Ruffini et al. 0). 


with very high accuracy. Figs. 12711301 shows the results for the luminosities in the 
30-400 keV and 2-10 keV bands, including the “prompt radiation”. Subsequently, 
the theoretically predicted GRB spectra have been evaluated at selected values of 
the arrival time (Ruffini et al. (2J)). 

The splendid news received the evening before the presentation of this talk 
is graphically represented by the XMM observations shown in Fig. 1201 Again, 
the XMM observations, like the corresponding ones of GRB 980425, occur after 
the decaying part of the afterglow and, in analogy to the one occurring in the 
system GRB 980425 / SN 1998bw / URCA-1, we call this source URCA-2. Further 
observations by XMM are highly recommended to follow the URCA-2 temporal 
evolution. Also in this system we are dealing with an astrophysical “triptych”: 
GRB 030329 / SN 2003dh / URCA-2. 
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Figure 27. The luminosity in the 2-10 keV and in the 30-400 keV energy bands predicted by our 
model are fitted to the data of R-XTE (GCN Circ. 1996 J98lR and HETE-2 (GCN Circ. 1997 
J99ID respectively. The SN 2003dh optical luminosity is given by the crosses (Hjorth et al. Eq 3)). 
Details in Bernardini et al. ©. 


10. Conclusions 

10.1. On the GRB-Supernova connection 

We first stress some general considerations originating from comparing and con¬ 
trasting the three GRB sources we have discussed: 

(1) The value of the B parameter for all three sources occurs, as theoretically 
expected, in the allowed range (see Fig. 1221) 

10“ 5 < B < 10~ 2 . (42) 

We have in fact: 

GRB 991216 B = 3.0 x 10“ 3 E dya = 4.8 x 10 53 erg 

GRB 980425 B = 7.0 x 10~ 3 E dya = 1.1 x 10 48 erg 

GRB 030329 B = 4.8 x 10~ 3 E dya = 2.1 x 10 52 erg 

(2) The enormous difference in the GRB energy of the sources simply relates 

to the electromagnetic energy of the black hole given in Eq. m which turns 
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Figure 28. The details of the theoretical fit of the prompt radiation of GRB 030329 have been 
reproduced by the filamentary structure in the ISM in our model. Details in Bernardini et al. 


out to be smaller than the critical value given by Eq. 0 . The fact that the 
theory is valid over 5 orders of magnitude is indeed very satisfactory. 

(3) Also revealing is the fact that in both sources GRB 980425 and GRB 030329 
the associated supernova energies are similar. We have, in fact, for both SN 
1998bw and SN 2003dh an energy ~ 10 49 erg. Details in Fraschetti et al. 
(lfi) and Bernardini et al. fQf>). The further comparison between the SN 
luminosity and the GRB intensity is crucial. In the case of GRB 980425 the 
GRB and the SN energies are comparable, and no dominance of one source 
over the other can be ascertained. In the case of GRB 030329 the energy of 
the GRB source is 10 3 larger than the SN: in no way the GRB can originate 
from the SN event. 

The above stringent energetics considerations and the fact that GRBs occur also 
without an observed supernova give a strong evidence that GRBs cannot originate 
from supernovae. 


10.2. URCA-1 and URCA-2 

We turn now to the most exciting search for the nature of URCA-1 and URCA-2. We 
have already mentioned above that a variety of possibilities naturally appear. The 
first possibility is that the URCA sources are related to the black hole originating 
the GRB phenomenon. In order to probe such an hypothesis, it would be very 
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Figure 29. The perfect fit of the late part of the afterglow of our theoretical model for the 2-10 
keV energy bands. The data refers to the R-XTE observations (GCN Circ. 1996 1 9 ti l ). Details in 
Bernardini et al. (96). 


important to find even a single case in which an URCA source occurs in association 
with a GRB and in absence of an associated supernova. Such a result, theoretically 
unexpected, would open an entire new problematic in relativistic astrophysics and 
in the physics of black holes. 

If indeed, as we expect, the clear association between URCA sources and the 
supernovae occurring together with the GRBs, then it is clear that the analysis 
of the other two possibilities will be favored. Namely, an emission from processes 
occurring in the early phases of the expansion of the supernova remnant or the very 
exciting possibility that for the first time we are observing a newly born neutron 
star out of the supernova phenomenon. Of course, this last hypothesis is the most 
important one, since it would offer new fundamental information about the outcome 
of the gravitational collapse, about the equations of state at supranuclear densities 
and about a variety of fundamental issues of relativistic astrophysics of neutron 
stars. We shall focus in the following only on this last topic. 

We have already recalled how the need for a rapid cooling process due to neu¬ 
trino anti-neutrino emission in the process of gravitational collapse leading to the 
formation of a neutron star was considered for the first time by George Gamow and 
Mario Schoenberg in 1941 B - It was Gamow who gave this process the name 
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GRB 030329 - SN2003dh - URCA-2 



Figure 30. The dotted line represents our theoretically predicted GRB030329 light curve in 7 - 
rays (30-400 keV) wit h the horizontal bar corresponding to the mean peak flux from HETE-2 
(GCN Circ. 1997 l99lB . The solid line represents the corresponding one in X-rays (2-10 keV) 
with the experimental data obtained by R-XTE (GCN Circ. 1996 (198m . The remaining points 
refer respectively to the op tical VLT data (Hjorth et al. 1 1 0(1 )) of SN2003bw and to the X-ray 
XMM data (Tiengo et al. 1 1D1 Vl of URCA-2. The dash-dotted lines corresponds to qualitative 
luminosity curves expected for URCA-2. It is interesting to compare and contrast these results 
with the ones for GRB980425/SN1998bw (see Fig. 3 in Ruffini et al. j 1 1) I. Details in Ruffini et 
al. Ja3). 


“Urea process”, sec|Appcndix A|and|Appcndix B| Since then, a systemat ic analysis 


103; 1 


Tsuruta 


of the theory o f neu tron star coolin g was advanced by Ts uruta 
and Cameron IToT Tsuruta et al. ( [l06 ) and by Canuto lll07tf . The coming of age 
of X-ray observatories such as Einstein (1978-1981), EXOSAT (1983-1986), ROSAT 
(1990-1998), and the contemporary missions of Chandra and XMM-Newton since 
1999 dramatically presented an observational situation establishing very embarrass¬ 
ing and stringent upper limits to the surface temperatu re of neutron stars in well 
known historical supernova remnants (see e.g. Romani ifioT . It was so that, for 
some remnants, notably SN 1006 and the Tycho supernova, the upper limits to 
the surface temperatures were significantly lower than the temperatures given by 


standard cooling times (see e.g. Romani ( j 108 11. Much of the theoretical works has 
been mainly directed, therefore, to find theoretical arguments in order to explain 
such low surface temperature T s ~ 0.5-1.0 x 10 6 K — embarrassingly low, when 
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compared to the initial hot 


plosion (see e.g. Romani (10 


10 11 K) birth of a neutron star in a supernova ex- 
). Some very important s teps in th is direction of 


research ha ve been r epresented by the works of Van Riper ll 1 Oil 1 1 1 (]) . Lattimer and 
his g roup JIB EH and by the most extensive work of Yakovlev and his group 
1111.31) . The youngest neutron star to be searched for using its thermal emission in 
this cont ext h as been the pulsar PSR J0205+6449 in 3C 58 (see e.g. Yakovlev and 
Pethick : l_K>i). which is 820 years old! Recently, evidence for the detection of ther¬ 
mal emission from the crab nebula pulsar was reported by Trumper (11141) which is, 
again, 951 years old. 

In the case of URCA-1 and URCA-2, we are exploring a totally different regime: 
the X-ray emission possibly from a recently born neutron star in the first days - 
months of its existence, where no observations have yet been performed and no 
embarrassing constraints upper limits on the surface temperature exist. The reason 
of approaching first the issue of the thermal emission from the neutron star surface 
is extremely important, since in principle it can give information on the equations 
of state in the core at supranuclear densities and on the detailed mechanism of 
the formation of the neutron star itself and the related neutrino emission. It is of 
course possible that the neutron star is initially fast rotating and its early emission 
is dominated by the magnetospheric emission or by accretion processes from the 
remnant which would overshadow the thermal emission. In that case a periodic 
signal related to the neutron star rotational period should in principle be observable 
in a close enough GRB source provided the suitable instrumentation from the Earth. 

The literature on young born neutron star is relatively scarse today. There are 
some very interesting contributions which state: “The time for a neutron star’s 
center to cool by the direct URCA process to a temperature T has been estimated 
to be t = 20 [T/ (10 9 AT)] 4 s. The direct URCA process and all the exotic cooling 
mechanisms only occur at supranuclear densities. Matter at subnuclear densities 
in neutron star crust cools primarily by diffusion of heat to the interior. Thus the 
surface temperature remains high, in the vicinity of 10 6 K or more, until the crust’s 
heat reservoir is consumed. After this diffusion time, which is on the order of 1 
100 years, the s urfac e temperature abruptly plunges to values below 5 x 10 5 K” 
(Lattimer et al. Il~2 )). “Soon after a supernova explosion, the young neutron star 
has large temperature gradients in the inner part of the crust. While the powerful 
neutrino emission quickly cools the core, the crust stays hot. The heat gradually 
flows inward on a conduction time scale and the whole process can be thought of 
as a cooling wave propagation from the center toward the surface” (Gnedin et al. 

EH))- 


The two considerations we have quoted above are developed in the case of spher¬ 
ical symmetry and we would like to keep the mind open, in this new astropliysical 
field, to additional factors, some more traditional than others, to be taken into 
account. Among the traditional ones we recall: 1) the presence of rotation and 
magnetic field which may affect the thermal conductivity and the structure of the 
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surface, as well as the above mentioned magnetosplieric emission; 2) there could 
be accretion of matter from the expanding nebula; and, among the nontraditional 
ones, we recall 3) some exciting theoretical possibilities advanced by Dyso n on vol¬ 
canoes on neutron stars EH as well as iron helide on neutron star JlT/t l. as well 
as the possib ility of piconuclear reactions on neutron star surface discussed in Lai 
& Salpeter Jlisf l. 

All the above are just scientific arguments to attract attention on the abrupt 
fall in luminosity reported in this meeting on URCA-1 by Elena Pian which is 
therefore, in this light, of the greatest scientific interest and further analysis should 
be followed to check if a similar behavior will be found in future XMM and Chandra 
observations also in URCA-2. 


10.3. Astrophysical implications 

In addition to these very rich problematics in the field of theoretical physics and 
theoretical astrophysics, there are also more classical astronomical and astrophysical 
issues, which will need to be answered if indeed the observations of a young neutron 
star will be confirmed. An important issue to be addressed will be how the young 
neutron star can be observed, escaping from being buried under the expelled matter 
of the collapsing star. A possible explanation can originate from the binary nature 
of the newly born neutron star: the binary system being formed by the newly 
formed black hole and the triggered gravitational collapse of a companion evolved 
star leading, possibly, to a “kick” on and ejection of the newly born neutron star. 
Another possibility, also related to the binary nature of the system, is that the 
supernova progenitor star has been depleted of its outer layer by dynamic tidal 
effects. 

In addition, there are other topics in which our scenario can open new research 
directions in fundamental physics and astrophysics: 

1) The problem of the instability leading to the complete gravitational collapse of a 
~ 1OM 0 star needs the introduction of a new critical mass for gravitational collapse, 
which is quite different from the one for white dwarfs and neutron stars which has 
been widely discussed in the current literature (see e.g. Giacconi & Ruffini (§)). 

2) The issue of the trigger of the instability of gravitational collapse induced by 
the GRB on the progenitor star of the supernova or, vice versa, by the supernova 
on the progenitor star of the GRB needs accurate timing and the considerations of 
new relativistic phenomena. 

3) The general relativistic instability induced on a nearby star by the formation of a 
black hole needs some very basic new developments in the field of general relativity. 

Only a very preliminary work exists on this subject, by Jim Wilson and his 
collaborators, see e.g. the paper by Mathews and Wilson in these proceedings 
ma>. The reason for the complexity in answering such a question is simply stated: 
unlike the majority of theoretical work on black holes, which deals mainly with one- 
body solutions, we have to address here a many-body problem in general relativity. 
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We are star ting in these days to recon sider, in this fram ework, some cl assic work by 
Fermi El, Hanni and Ruffini ill 2 il l. Majumdar Il22l i. Papapetrou (1 12. ‘ill . Parker 
et al. (124), Bird et al. ( 1125 1 which may lead to a new understanding of general 
relativistic effects relevant to these astrophysical “triptychs”. 
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Appendix A. On the Urea process 

From G. Gamow m 

“The summer of 1939 I spent with my family vacation on the Copacabana beach 
in Rio de Janeiro. One evening, visiting the famous Casino da Urea to watch the 
gamblers, I was introduced to a young theoretical physicist born on an Amazon 
River plantation, named Mario Schoenberg. We became friends, and I arranged for 
him a Guggenheim fellowship to spend a year in Washington to work with me in 
nuclear astrophysics. His visit was very successful, and we hit upon a process which 
could be responsible for the vast stellar explosions known as supernovae. The trick 
is done by alternative absorption and reemission of one of the thermal electrons in 
the very hot (billions of degrees!) stellar interior by various atomic nucleai. Both 
processes are accompanied by the emission of neutrinos and antineutrinos which, 
possessing tremendous penetrating power, pass through the body of a star like a 
swarm of musquitoes through chicken wire and carry with them large amount of 
energy. Thus, the stellar interior cools rapidly, the pressure drops, and the stellar 
body collapse with a great explosion of light and heat. 

All this is too complicated to explain in nontechnical words, and I am mentioning 
it only as background for how we came to give that process its name. We called 
it the Urea process, partially to commemorate the casino in which we first met, 
and partially because the Urea process results in a rapid disappearance of thermal 
energy from the interior of the star, similar to the rapid disappearance of money 
from the pockets of the gamblers of the Casino da Urea. Sending our article “On 
the Urea process” for publication in The Physical Review I was worried that the 
Editor would ask why we called the process “Urea”. After much thought I decided 
to say that this is short for “UnRecordable Cooling Agent”, but they never asked. 
Today, there are other known cooling processes involving neutrinos which work even 
faster than the Urea process. For example, a neutrino pair can be formed instead 
of two gamma quanta in the annihilation of a positive and negative electrons”. 
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